Genomic epidemiology reveals geographical clustering of multidrug-resistant Escherichia coli ST131 associated with bacteraemia in Wales.

Antibiotic resistance is a significant global public health concern. Uropathogenic Escherichia coli sequence type (ST)131, a widely prevalent multidrug-resistant clone, is frequently associated with bacteraemia. This study investigates third-generation cephalosporin resistance in bloodstream infections caused by E. coli ST131. From 2013-2014 blood culture surveillance in Wales, 142 E. coli ST131 genomes were studied alongside global data. All three major ST131 clades were represented across Wales, with clade C/H30 predominant (n = 102/142, 71.8%). Consistent with global findings, Welsh strains of clade C/H30 contain ß-lactamase genes from the blaCTX-M-1 group (n = 65/102, 63.7%), which confer resistance to third-generation cephalosporins. Most Welsh clade C/H30 genomes belonged to sub-clade C2/H30Rx (58.3%). A Wales-specific sub-lineage, named GB-WLS.C2, diverged around 1996-2000. An introduction to North Wales around 2002 led to a localised cluster by 2009, depicting limited genomic diversity within North Wales. This investigation emphasises the value of genomic epidemiology, allowing the detection of genetically similar strains in local areas, enabling targeted and timely public health interventions.

A Proofreading Mutation with an Allosteric Effect Allows a Cluster of SARS-CoV-2 Viruses to Rapidly Evolve.

The RNA-dependent RNA polymerase of the severe acute respiratory syndrome coronavirus 2 virus is error prone, with errors being corrected by the exonuclease (NSP14) proofreading mechanism. However, the mutagenesis and subsequent evolutionary trajectory of the virus is mediated by the delicate interplay of replicase fidelity and environmental pressures. Here, we have shown that a single, distal mutation (F60S) in NSP14 can have a profound impact upon proofreading with an increased accumulation of mutations and elevated evolutionary rate being observed. Understanding the implications of these changes is crucial, as these underlying mutational processes may have important implications for understanding the population-wide evolution of the virus. This study underscores the urgent need for continued research into the replicative mechanisms of this virus to combat its continued impact on global health, through the re-emergence of immuno-evasive variants.

Identification of two distinct phylogenomic lineages and model strains for the understudied cystic fibrosis lung pathogen Burkholderia multivorans.

Burkholderia multivorans is the dominant Burkholderia pathogen recovered from lung infection in people with cystic fibrosis. However, as an understudied pathogen there are knowledge gaps in relation to its population biology, phenotypic traits and useful model strains. A phylogenomic study of B. multivorans was undertaken using a total of 283 genomes, of which 73 were sequenced and 49 phenotypically characterized as part of this study. Average nucleotide identity analysis (ANI) and phylogenetic alignment of core genes demonstrated that the B. multivorans population separated into two distinct evolutionary clades, defined as lineage 1 (n=58 genomes) and lineage 2 (n=221 genomes). To examine the population biology of B. multivorans, a representative subgroup of 77 B. multivorans genomes (28 from the reference databases and the 49 novel short-read genome sequences) were selected based on multilocus sequence typing (MLST), isolation source and phylogenetic placement criteria. Comparative genomics was used to identify B. multivorans lineage-specific genes - ghrB_1 in lineage 1 and glnM_2 in lineage 2 - and diagnostic PCRs targeting them were successfully developed. Phenotypic analysis of 49 representative B. multivorans strains showed considerable inter-strain variance, but the majority of the isolates tested were motile and capable of biofilm formation. A striking absence of B. multivorans protease activity in vitro was observed, but no lineage-specific phenotypic differences were demonstrated. Using phylogenomic and phenotypic criteria, three model B. multivorans CF strains were identified, BCC0084 (lineage 1), BCC1272 (lineage 2a) and BCC0033 lineage 2b, and their complete genome sequences determined. B. multivorans CF strains BCC0033 and BCC0084, and the environmental reference strain, ATCC 17616, were all capable of short-term survival within a murine lung infection model. By mapping the population biology, identifying lineage-specific PCRs and model strains, we provide much needed baseline resources for future studies of B. multivorans.

Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Novel recombinant SARS-CoV-2 lineage detected through genomic surveillance in Wales, UK.

Recombination, the process whereby a segment of genetic material from one genome is inserted into another, producing a new chimeric genome, is an important evolutionary mechanism frequently observed in coronaviruses. The risks posed by recombination include the shuffling of advantageous mutations that may increase transmissibility, severity or vaccine escape. We present a genomic and epidemiological description of a new recombinant lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), XR, first identified in Wales. The Pathogen Genomics Unit (Public Health Wales, UK) sequences positive SARS-CoV-2 PCR tests using the ARTIC SARS-CoV-2 sequencing protocol. Recombinants were detected using an in-house pipeline and the epidemiological data analysed in R. Nosocomial cases were defined as those with samples taken after >7 days in hospital. Between February and March 2022, we identified 78 samples with highly similar genomes, comprising a BA.1-like 5' end, a BA.2-like 3' end and a BA.2-like spike protein. This signature is consistent with recombination and was defined as XR by Pangolin (PANGO v1.8). A total of 50 % of cases had a sample collected whilst in hospital and the first three cases were immunocompromised patients. The patient median age was 58 years (range: 4-95 years) and most of the patients were fully vaccinated against SARS-CoV-2 (74 % third dose/booster). Three patients died within 28 days of their sample collection date, one of whom had COVID-19 listed amongst ICD10 (International Classification of Diseases 10) coded causes of death. Our integrated system enabled real-time monitoring of recombinant SARS-CoV-2 for early detection, in order to rapidly risk assess and respond. This work highlights the importance of setting-based surveillance of recombinant SARS-CoV-2, as well as the need to monitor immunocompromised populations through repeat testing and sequencing.

Emergence of immune escape at dominant SARS-CoV-2 killer T cell epitope.

We studied the prevalent cytotoxic CD8 T cell response mounted against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Spike glycoprotein269-277 epitope (sequence YLQPRTFLL) via the most frequent human leukocyte antigen (HLA) class I worldwide, HLA A∗02. The Spike P272L mutation that has arisen in at least 112 different SARS-CoV-2 lineages to date, including in lineages classified as "variants of concern," was not recognized by the large CD8 T cell response seen across cohorts of HLA A∗02+ convalescent patients and individuals vaccinated against SARS-CoV-2, despite these responses comprising of over 175 different individual T cell receptors. Viral escape at prevalent T cell epitopes restricted by high frequency HLAs may be particularly problematic when vaccine immunity is focused on a single protein such as SARS-CoV-2 Spike, providing a strong argument for inclusion of multiple viral proteins in next generation vaccines and highlighting the need for monitoring T cell escape in new SARS-CoV-2 variants.

Epidemiological analysis of the first 1000 cases of SARS-CoV-2 lineage BA.1 (B.1.1.529, Omicron) compared with co-circulating Delta in Wales, UK.

BACKGROUND: The Omicron (lineage B.1.1.529) variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported in Wales, UK, on 3 December 2021. The aim of the study was to describe the first 1000 cases of the Omicron variant by demographic, vaccination status, travel and severe outcome status and compare this to contemporaneous cases of the Delta variant. METHODS: Testing, typing and contact tracing data were collected by Public Health Wales and analysis undertaken by the Communicable Disease Surveillance Centre (CDSC). Risk ratios for demographic factors and symptoms were calculated comparing Omicron cases to Delta cases identified over the same time period. RESULTS: By 14 December 2021, 1000 cases of the Omicron variant had been identified in Wales. Of the first 1000, just 3% of cases had a prior history of travel revealing rapid community transmission. A higher proportion of Omicron cases were identified in individuals aged 20-39, and most cases were double vaccinated (65.9%) or boosted (15.7%). Age-adjusted analysis also revealed that Omicron cases were less likely to be hospitalised (0.4%) or report symptoms (60.8%). Specifically a significant reduction was observed in the proportion of Omicron cases reporting anosmia (8.9%). CONCLUSION: Key findings include a lower risk of anosmia and a reduced risk of hospitalisation in the first 1000 Omicron cases compared with co-circulating Delta cases. We also identify that existing measures for travel restrictions to control importations of new variants identified outside the United Kingdom did not prevent the rapid ingress of Omicron within Wales.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic.

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Future-proofing and maximizing the utility of metadata: The PHA4GE SARS-CoV-2 contextual data specification package.

BACKGROUND: The Public Health Alliance for Genomic Epidemiology (PHA4GE) (https://pha4ge.org) is a global coalition that is actively working to establish consensus standards, document and share best practices, improve the availability of critical bioinformatics tools and resources, and advocate for greater openness, interoperability, accessibility, and reproducibility in public health microbial bioinformatics. In the face of the current pandemic, PHA4GE has identified a need for a fit-for-purpose, open-source SARS-CoV-2 contextual data standard. RESULTS: As such, we have developed a SARS-CoV-2 contextual data specification package based on harmonizable, publicly available community standards. The specification can be implemented via a collection template, as well as an array of protocols and tools to support both the harmonization and submission of sequence data and contextual information to public biorepositories. CONCLUSIONS: Well-structured, rich contextual data add value, promote reuse, and enable aggregation and integration of disparate datasets. Adoption of the proposed standard and practices will better enable interoperability between datasets and systems, improve the consistency and utility of generated data, and ultimately facilitate novel insights and discoveries in SARS-CoV-2 and COVID-19. The package is now supported by the NCBI's BioSample database.

Novel application of metagenomics for the strain-level detection of bacterial contaminants within non-sterile industrial products - a retrospective, real-time analysis.

The home and personal care (HPC) industry generally relies on initial cultivation and subsequent biochemical testing for the identification of microorganisms in contaminated products. This process is slow (several days for growth), labour intensive, and misses organisms which fail to revive from the harsh environment of preserved consumer products. Since manufacturing within the HPC industry is high-throughput, the process of identification of microbial contamination could benefit from the multiple cultivation-independent methodologies that have developed for the detection and analysis of microbes. We describe a novel workflow starting with automated DNA extraction directly from a HPC product, and subsequently applying metagenomic methodologies for species and strain-level identification of bacteria. The workflow was validated by application to a historic microbial contamination of a general-purpose cleaner (GPC). A single strain of Pseudomonas oleovorans was detected metagenomically within the product. The metagenome mirrored that of a contaminant isolated in parallel by a traditional cultivation-based approach. Using a dilution series of the incident sample, we also provide evidence to show that the workflow enables detection of contaminant organisms down to 100 CFU/ml of product. To our knowledge, this is the first validated example of metagenomics analysis providing confirmatory evidence of a traditionally isolated contaminant organism, in a HPC product.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic.

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

Discovery of the Pseudomonas Polyyne Protegencin by a Phylogeny-Guided Study of Polyyne Biosynthetic Gene Cluster Diversity.

Natural products that possess alkyne or polyyne moieties have been isolated from a variety of biological sources and possess a broad a range of bioactivities. In bacteria, the basic biosynthesis of polyynes is known, but their biosynthetic gene cluster (BGC) distribution and evolutionary relationship to alkyne biosynthesis have not been addressed. Through comprehensive genomic and phylogenetic analyses, the distribution of alkyne biosynthesis gene cassettes throughout bacteria was explored, revealing evidence of multiple horizontal gene transfer events. After investigation of the evolutionary connection between alkyne and polyyne biosynthesis, a monophyletic clade was identified that possessed a conserved seven-gene cassette for polyyne biosynthesis that built upon the conserved three-gene cassette for alkyne biosynthesis. Further diversity mapping of the conserved polyyne gene cassette revealed a phylogenetic subclade for an uncharacterized polyyne BGC present in several Pseudomonas species, designated pgn. Pathway mutagenesis and high-resolution analytical chemistry showed the Pseudomonas protegens pgn BGC directed the biosynthesis of a novel polyyne, protegencin. Exploration of the biosynthetic logic behind polyyne production, through BGC mutagenesis and analytical chemistry, highlighted the essentiality of a triad of desaturase proteins and a thioesterase in both the P. protegens pgn and Trinickia caryophylli (formerly Burkholderia caryophylli) caryoynencin pathways. We have unified and expanded knowledge of polyyne diversity and uniquely demonstrated that alkyne and polyyne biosynthetic gene clusters are evolutionarily related and widely distributed within bacteria. The systematic mapping of conserved biosynthetic genes across the available bacterial genomic diversity proved to be a fruitful method for discovering new natural products and better understanding polyyne biosynthesis. IMPORTANCE Natural products bearing alkyne (triple carbon bond) or polyyne (multiple alternating single and triple carbon bonds) moieties exhibit a broad range of important biological activities. Polyyne metabolites have been implicated in important ecological roles such as cepacin mediating biological control of plant pathogens and caryoynencin protecting Lagriinae beetle eggs against pathogenic fungi. After further phylogenetic exploration of polyyne diversity, we identified a novel gene cluster in Pseudomonas bacteria with known biological control abilities and proved it was responsible for synthesizing a new polyyne metabolite, protegencin. The evolutionary analysis of polyyne pathways showed that multiple biosynthetic genes were conserved, and using mutagenesis, their essentiality was demonstrated. Our research provides a foundation for the future modification of polyyne metabolites and has identified a novel polyyne, protegencin, with potential bioactive roles of ecological and agricultural importance.

CLIMB-COVID: continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance.

In response to the ongoing SARS-CoV-2 pandemic in the UK, the COVID-19 Genomics UK (COG-UK) consortium was formed to rapidly sequence SARS-CoV-2 genomes as part of a national-scale genomic surveillance strategy. The network consists of universities, academic institutes, regional sequencing centres and the four UK Public Health Agencies. We describe the development and deployment of CLIMB-COVID, an encompassing digital infrastructure to address the challenge of collecting and integrating both genomic sequencing data and sample-associated metadata produced across the COG-UK network.

Monitoring SARS-CoV-2 in municipal wastewater to evaluate the success of lockdown measures for controlling COVID-19 in the UK.

SARS-CoV-2 and the resulting COVID-19 pandemic represents one of the greatest recent threats to human health, wellbeing and economic growth. Wastewater-based epidemiology (WBE) of human viruses can be a useful tool for population-scale monitoring of SARS-CoV-2 prevalence and epidemiology to help prevent further spread of the disease, particularly within urban centres. Here, we present a longitudinal analysis (March-July 2020) of SARS-CoV-2 RNA prevalence in sewage across six major urban centres in the UK (total population equivalent 3 million) by q(RT-)PCR and viral genome sequencing. Our results demonstrate that levels of SARS-CoV-2 RNA generally correlated with the abundance of clinical cases recorded within the community in large urban centres, with a marked decline in SARS-CoV-2 RNA abundance following the implementation of lockdown measures. The strength of this association was weaker in areas with lower confirmed COVID-19 case numbers. Further, sequence analysis of SARS-CoV-2 from wastewater suggested that multiple genetically distinct clusters were co-circulating in the local populations covered by our sample sites, and that the genetic variants observed in wastewater reflected similar SNPs observed in contemporaneous samples from cases tested in clinical diagnostic laboratories. We demonstrate how WBE can be used for both community-level detection and tracking of SARS-CoV-2 and other virus' prevalence, and can inform public health policy decisions. Although, greater understanding of the factors that affect SARS-CoV-2 RNA concentration in wastewater are needed for the full integration of WBE data into outbreak surveillance. In conclusion, our results lend support to the use of routine WBE for monitoring of SARS-CoV-2 and other human pathogenic viruses circulating in the population and assessment of the effectiveness of disease control measures.

Genomics reveals the novel species placement of industrial contaminant isolates incorrectly identified as Burkholderia lata.

The Burkholderia cepacia complex (Bcc) is a closely related group of bacteria, composed of at least 20 different species, the accurate identification of which is essential in the context of infectious diseases. In industry, they can contaminate non-food products, including home and personal care products and cosmetics. The Bcc are problematic contaminants due to their ubiquitous presence and intrinsic antimicrobial resistance, which enables them to occasionally overcome preservation systems in non-sterile products. Burkholderia lata and Burkholderia contaminans are amongst the Bcc bacteria encountered most frequently as industrial contaminants, but their identification is not straightforward. Both species were historically established as a part of a group known collectively as taxon K, based upon analysis of the recA gene and multilocus sequence typing (MLST). Here, we deploy a straightforward genomics-based workflow for accurate Bcc classification using average nucleotide identity (ANI) and core-gene analysis. The workflow was used to examine a panel of 23 Burkholderia taxon K industrial strains, which, based on MLST, comprised 13 B. lata, 4 B. contaminans and 6 unclassified Bcc strains. Our genomic identification showed that the B. contaminans strains retained their classification, whilst the remaining strains were reclassified as Burkholderia aenigmatica sp. nov. Incorrect taxonomic identification of industrial contaminants is a problematic issue. Application and testing of our genomic workflow allowed the correct classification of 23 Bcc industrial strains, and also indicated that B. aenigmatica sp. nov. may have greater importance than B. lata as a contaminant species. Our study illustrates how the non-food manufacturing industry can harness whole-genome sequencing to better understand antimicrobial-resistant bacteria affecting their products.

SARS-CoV-2 within-host diversity and transmission.

Extensive global sampling and sequencing of the pandemic virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have enabled researchers to monitor its spread and to identify concerning new variants. Two important determinants of variant spread are how frequently they arise within individuals and how likely they are to be transmitted. To characterize within-host diversity and transmission, we deep-sequenced 1313 clinical samples from the United Kingdom. SARS-CoV-2 infections are characterized by low levels of within-host diversity when viral loads are high and by a narrow bottleneck at transmission. Most variants are either lost or occasionally fixed at the point of transmission, with minimal persistence of shared diversity, patterns that are readily observable on the phylogenetic tree. Our results suggest that transmission-enhancing and/or immune-escape SARS-CoV-2 variants are likely to arise infrequently but could spread rapidly if successfully transmitted.

Kill and cure: genomic phylogeny and bioactivity of Burkholderia gladioli bacteria capable of pathogenic and beneficial lifestyles.

Burkholderia gladioli is a bacterium with a broad ecology spanning disease in humans, animals and plants, but also encompassing multiple beneficial interactions. It is a plant pathogen, a toxin-producing food-poisoning agent, and causes lung infections in people with cystic fibrosis (CF). Contrasting beneficial traits include antifungal production exploited by insects to protect their eggs, plant protective abilities and antibiotic biosynthesis. We explored the genomic diversity and specialized metabolic potential of 206 B. gladioli strains, phylogenomically defining 5 clades. Historical disease pathovars (pv.) B. gladioli pv. allicola and B. gladioli pv. cocovenenans were distinct, while B. gladioli pv. gladioli and B. gladioli pv. agaricicola were indistinguishable; soft-rot disease and CF infection were conserved across all pathovars. Biosynthetic gene clusters (BGCs) for toxoflavin, caryoynencin and enacyloxin were dispersed across B. gladioli, but bongkrekic acid and gladiolin production were clade-specific. Strikingly, 13 % of CF infection strains characterized were bongkrekic acid-positive, uniquely linking this food-poisoning toxin to this aspect of B. gladioli disease. Mapping the population biology and metabolite production of B. gladioli has shed light on its diverse ecology, and by demonstrating that the antibiotic trimethoprim suppresses bongkrekic acid production, a potential therapeutic strategy to minimize poisoning risk in CF has been identified.

Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK.

The United Kingdom's COVID-19 epidemic during early 2020 was one of world's largest and was unusually well represented by virus genomic sampling. We determined the fine-scale genetic lineage structure of this epidemic through analysis of 50,887 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes, including 26,181 from the UK sampled throughout the country's first wave of infection. Using large-scale phylogenetic analyses combined with epidemiological and travel data, we quantified the size, spatiotemporal origins, and persistence of genetically distinct UK transmission lineages. Rapid fluctuations in virus importation rates resulted in >1000 lineages; those introduced prior to national lockdown tended to be larger and more dispersed. Lineage importation and regional lineage diversity declined after lockdown, whereas lineage elimination was size-dependent. We discuss the implications of our genetic perspective on transmission dynamics for COVID-19 epidemiology and control.

A structure-based nomenclature for Bacillus thuringiensis and other bacteria-derived pesticidal proteins.

In 1998 a nomenclature for the growing list of pesticidal proteins from Bacillus thuringiensis (Bt) was derived based solely on protein sequence comparisons. This nomenclature was widely adopted and provided a robust framework for the naming and classification of the proteins. The success of these proteins in integrated pest management schemes prompted an increased effort to find others with improved or more diverse activities. These discovery activities led to the characterization of proteins from a wider range of bacteria and with a variety of different protein folds. Since most of these new proteins were grouped together as Cry proteins it became apparent that the existing nomenclature had limitations in representing the diverse range of proteins that had been identified. This revised nomenclature retains the basic principles of the 1998 version but provides specific mnemonics to represent different structural groups. For the purposes of consistency, the vast majority of the proteins have either retained their name or have a new name that clearly references the previous one. Other pesticidal proteins not previously included in the nomenclature have been incorporated into this version.

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.